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Abstract 

This paper deals with the challenging problem of spectrum sensing in cognitive radio. We consider a 
stochastic system model where the Primary User (PU) transmits a periodic signal over fading channels. 
The effect of frequency offsets due to oscillator mismatch, and Doppler offset is studied. We show 
that for this case the Likelihood Ratio Test (LRT) cannot be evaluated pointwise. We present a novel 
approach to approximate the marginilisation of the frequency offset using a single point estimate. This 
is obtained via a low complexity Constrained Adaptive Notch Filter (CANF) to estimate the frequency 
offset. Performance is evaluated via numerical simulations and it is shown that the proposed spectrum 
sensing scheme can achieve the same performance as the "near-optimal" scheme, that is based on a 
bank of matched filters, using only a fraction of the complexity required. 

Index Terms 

Spectrum Sensing, Likelihood Ratio Test, Adaptive Notch Filter. 

I. Introduction 

In recent years, cognitive radio has attracted intensive research focus due to the pressing 
demand for efficient usage of the frequency spectrum [IJ. In a cognitive radio system, secondary 
radio users try to find "blank spaces", in which the licensed frequency band is not being used by 
primary radio users, for communications. A key problem in cognitive radio is that the secondary 
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users need to vacate the frequency band as quickly as possible if the corresponding primary 
radio emerges, and begins transmission. 

Spectrum sensing is a mandatory functionality in any CR-based wireless system that shares 
spectrum bands with primary services, such as the IEEE 802.22 standard which proposes 
to reuse vacant spectrum in the TV broadcast bands. There has been significant research on 
spectrum sensing for cognitive radio, see |l3l, JH for an overview. Essentially, spectrum sensing 
is a decision making or classification problem, of the form involving first an estimation stage, 
followed a decision stage. The secondary network needs to make a decision between two possible 
hypotheses given an observation vector: that the frequency band is either occupied or vacant. 
The more knowledge we have on the nature of the primary user's signal, the more reliable our 
decision process will become. If no prior knowledge of the primary user's signal, the energy 
detector based approaches (also called radiometry) are the most common for spectrum sensing, 
mainly due to their low computational complexity, see [5J and references within. If additional 
prior knowledge about the primary signal exists, more advanced techniques can be explored. For 
example, if the primary signal is a priori known, the matched filter is optimal in the sense that 
it maximises the Signal to Noise Ratio (SNR) [[51, jU. 

In practical systems, such a priori knowledge often exists in the form of a pilot signal, which 
is used by the primary network, enabling a Waveform-Based Sensing For example, in the 
standard of digital TV (DTV) system, made by advanced television standard committee (ATSC), 
there are multiple sinusoid pilots located at different frequencies 

Previous works: A few papers considered the problem of spectrum sensing using a pilot 
signal in the form of sinusoids. These include |I3, |[8l and ^ where the detection scheme was 
designed for the case that the primary user transmits at a known deterministic frequency and 
the amplitude/channel may or may not be known a priori. 

Contribution: in contrast to those papers, we extend the system model and consider two practical 
effects which are of significance in a wireless communication system: first, we consider the 
case of unknown Rayleigh fading channel. Second, we allow for frequency offsets due to 
Doppler offset and mismatched oscillators being present in the communications system. These 
two assumptions make the sensing algorithm complex. In order to overcome this difficulty, we 
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shall present two "close to optimal" sensing algorithms that are based on a bank of matched 
filters, and on a frequency estimation (periodogram) approaches. While these solutions perform 
close to optimal, their high complexity may prevent them from being practical. Instead, we 
propose a low complexity algorithm that is based on the principle of Generalised Likelihood 
Ratio Test (GLRT), where we condition the LRT on an estimate of the nuisance frequency 
parameter (the estimate is not the Maximum Likelihood (ML) estimate, and therefore, the 
proposed algorithm is not the GLRT). The estimation of the nuisance frequency parameter is 
based on the Constrained Adaptive Notch Filter (CANF). Then, we perform a single matched 
filter centered on the estimated frequency, providing very good performance with only a fraction 
of the computational complexity. 

The following notation is used throughout: random variables are denoted by upper case letters 
and their realizations by lower case letters. In addition, bold will be used to denote a vector or 
matrix quantity. 

II. System Model 

Consider a secondary radio communication network using a single sensor. The secondary 
sensor needs to monitor the activity of the primary network. We assume that the primary radio 
signal uses a pilot signal in the form of a sinusoid, as in [|71, |[8l and lIH, see Fig[T] For example, 
in High-definition TV (HDTV) systems, developed by advanced television standard committee 
(ATSC), there are multiple sinusoid pilots located at different frequency points. 



A. Model Assumptions 

We introduce the system model for spectrum sensing 

• Assume a secondary radio communication network using a single receive antenna. 

• The primary is active in a frame (block) of length M samples, with probability P (Ho) and 
idle with probability P{Ui). 

• The observation at the sensor can be written as the following binary hypothesis: 
no ■■ Y{m) = V{m), m = l,...,M 
Ui : Y{m) = H sin {niQ + 0) + V{m), m = 1, . . . , M, 
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where H is the unknown random channel gain (assumed to be fixed throughout a frame), 
Vt is the fixed unknown random angular frequency of the primary user signal, which is 
not known a-priori, due to oscillator mismatch and Doppler offset, and is the random 
unknown phase offset of the pilot signal, and V{m) is the random additive noise. 

B. Prior specification 

We present the relevant aspects of the Bayesian model. 

• The channel gain H is Rayleigh fading. Therefore its density function is 



where cr^ is the known variance of the channel gain. 

• The phase offset O is random and follows a uniform density function O ~ [/ [0, 27r) . 

• The angular frequency of the primary user's signal, Q, can be parametrised by density 
function distribution, p (u), over the support \uJ — e,uJ + e], with ZJ being the nominal 
known angular frequency and e as the maximal offset, determined by the Doppler offset and 
oscillators mismatch. In this paper, for ease of presentation, we assume a uniform prior 

• The received signal is corrupt by zero-mean i.i.d. additive white Gaussian noise (AWGN) 
V{m) ~ CN {0, (7^), with a known variance and power spectral density (PSD) Nq. 



The objective of spectrum sensing is to make a decision on the binary hypothesis testing 
(choose "Ho or "Hi) based on the received signal. 

A. Decision criterion 

Here we formulate the problem using Bayes' criterion flU. In doing so, two assumptions are 
made. First, the probabilities P ("Hq) and P (J-Li) are known. The second assumption is that a 
cost, Cxy, is assigned to each possible decision. Cxy is the associated cost of making a decision 
Hx, given that the true hypothesis is Hy. 




(2) 



in. Spectrum Sensing - Problem Definition 
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C = P{'Ho)[Coo p(yi:A/|'Ho) C?yi:A/ + Cio / ^'(yi:A/|'Ho)(iyi:A/ 

+ -P('Hi) I Coi / p{yi:M\'Hi)dyi:M + Cn p{yi:M\Hi) dyi;M 

\ JAo JAi 



The problem of designing the decision rule can be modeled as an optimization problem whose 
objective is to minimize the cost function 



(3) 



It can be shown that the optimum decision rule is a likelihood-ratio test given by 

■Hi 

A /V X A p(yi:A/l'^i) ^ P {^^o) Cio — Coo a ... 

^ ^^'-''^ - Ju;:i;\no) ^< pmco,-c,, - ^' 

where Cxy is the associated cost of making a decision T-L^, given that the true hypothesis is T-Ly, 
and we define yi.,M = [y (1) , . . . , y (Af)]. 

The major difficulty in using the LRT is its requirement on obtaining the exact distributions 
under each hypothesis in ([T]). Under the NULL hypothesis, the observations are independent, 
and the evaluation of the evidence, p (yi.^/l'Ho)^ can be decomposed as 

M 1 / 1 \ 

p{yi:M\no) = np(?/Mi^o) = -J— n ( -^^'(^) ) • 

m=l V V rn=l \ v / 

The distribution of the alternative, p(yi:A/|^i), may be harder to obtain depending on the 
knowledge of the system parameters. Here we develop the solution for several cases with different 
levels of knowledge of the system parameters. 

B. Case I: no frequency offset (u = uj) and known channel gain H 

With no frequency offset present (tu is known exactly) and the channel gain h is known 

a-priori, the likelihood ratio has the following expression [[6l 

, p(yi:A/|Hi) [ -Mh^ \ (2h \ 

A Yi:A/ = —, r— = exp ^ Mo — r , (6) 

l>(yi:A/|Ho) \ 2al J \al J 

where Jq is the modified Bessel function and r is defined as 

r (cj) = ^/(yf+yf), where (7a) 

M 

Uc (w) = ^ y{m) cos {mu) , (7b) 

m=l 
M 

Us (w) = yim) sin (mw) . (7c) 

m.=l 
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Here a single matched filter is required to perform the LRT. 

C. Case II: no frequency offset (u = uJ) and unknown channel gain H 

With no frequency offset present, but the realisation of the channel gain h unknown a-priori, 
the LRT can be evaluated exactly by marginalising Q over the unknown channel,//, as 

j'OO 

A"(Yi.m)= / K'{Y,.,M\h)p{h)dh 
Jo 



al ( la, 



^^H^J'^U^J^^^H"^''' 



where r (cu) is defined in (TTal) . Again, a single matched filter is required to perform the LRT. 

D. Case III: unknown frequency and unknown channel gain 

When a frequency offset is present and the channel gain h is unknown a-priori, we consider 
([8]) and marginalise over the unknown random frequency u, 

(9) 



The integral in ([91) is not analytic and requires approximation techniques, and we shall derive 
numerical approximations of ^ in the following sections. 

IV. High-Complexity Blind Spectrum Sensing Estimation 

In this section we briefly present two possible approximations which have high complexity 
relative to the solution we propose in this paper. They act as benchmarks for comparison with 
our solution. The first approach is based on a bank of matched filters. We approximate p {u) 
using a discrete density function pd (w) with K discrete values as 

K 

Pd{uj) ^^p{ujk) S - ^k) , (10) 

k=l 
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where Uk = (uj - e) + kA, k = {1,2, K} , K = ^. 

We can now approximate dH) using a discretization, written as 

w ' a \ w \ w ' a) / 

(11) 

Here, the LR function is evaluated by using a bank of K matched filters, one per frequency w^. 
This solution has high complexity since it performs a weighted average of Clearly, with this 
approach a strong trade-off between performance and complexity burden occurs . The matched 
filter is sensitive to frequency mismatch (see Section|VI]for analysis) and it is therefore desirable 
to set K to be very large. This would make the gap between consecutive discrete frequencies, 
Acj, very small and make the frequency mismatch vs. w^) negligible. This however would 
result in a very costly implementation. In cases where only a few matched filters are used, it is 
likely that a frequency mismatch will occur, leading to poor performance. 

The second approach is based on a GLRT and discrete Fourier transform. This solution 
first produces a periodogram, followed by a LRT conditioned on the maximum value of the 
frequency obtained from the periodogram. This has high complexity due to the construction of 
the periodogram. In the GLRT, we condition on frequency 1) obtained from a point estimate 
from the periodogram [|6l. Hence, we approximate dH) as 



A'^' (Yi.m) ^ = , , exp " kl^A . (12) 

w ' a J / 

It is well known that Qml can be asymptotically obtained by maximising the periodogram H, so 



that fluL = argmax^ 



Em=i y M exp" 



'jLum 



where j = V — 1. The accuracy of the frequency 
estimator depends on the number of samples in the frame, M. As with the bank of matched 
filters, a fine grid of frequencies is required, resulting in a highly computational algorithm. If, 
on the other hand, we used only a coarse grid with only a few frequencies, that would result in 
a significant loss of accuracy and high estimation error. 

V. Novel Low Complexity Blind Spectrum Sensing Estimation 

In this section we present a novel algorithm to perform a low complexity spectrum sensing 
for Case III. As in the case of approach two, our solution is also based on the GLRT, but 
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replaces the grid search required by the periodogram construction with an adaptive notch fiker 
based frequency estimator. This reduces the computational complexity significantly for the same 
estimation accuracy. Here we develop a non-standard solution for the GLRT which involves 
designing a notch filter which performs adaptive frequency estimation. The key to our solution 
is to utilise the result involving the representation of the transmitted signal under "Hi in ([T]) as 
a 2-nd order autoregressive process, obtained via a trigonometric identity, 

sin {mfl + 6) = 2 cos (Q) sin ((m - l)n + 9) + sin ((m - 2)fi + 6) . (13) 

In the frequency domain, the transmitted signal is represented by Dirac masses, with unknown 
locations. Using (fT3l) . we can estimate this location via a localised filter in the family of notch 
filters. A notch filter is a filter that contains a null in its frequency response characteristics. Here, 
for simplicity, we concentrate on a 2-nd order Infinite Impulse Response (IIR) which contains 
a pair of complex-conjugate zeros on the unit circle and a pair of complex-conjugate poles at 
the same frequency inside the unit circle, and has the following transfer function 

1 — p{m)P{m)z~^ + p'^{m)z~'^ 
where the values /3(m) determine the centre of the notch filter frequency, and < p(m) < 1 
defines the location of the poles inside the unit disk. This design has the properties of having a 
symmetric frequency response and a narrow bandwidth, provided that p(m) is close to 1. This 
filter is simple to design, requiring the estimation of two parameters, keeping complexity low, 
whilst providing narrow-band frequency selectivity. We now present the specific details of the 
proposed frequency estimation algorithm. 

The output of the filter, s(m), as defined by (fT4l) can be expressed as 

s(m) = y(m) + /3{m)y{m — 1) + y{m — 2) — p{m)(3{m)s{m — 1) — p'^{m)s{m — 2). (15) 

We formulate the joint optimisation for (3 and p using the following criterion 

1 / 1 \ 

B ■ < B<B '"=^ (16) 

P < Pmax 
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where /3min, /^max are determined by the allowed frequency offset e, and pmax < 1 to ensure the 

filter is stable, according to Lyapunov stability criterion IfTOll . This cost function is designed to 

achieve two goals: the first is to minimise the variance at the output of the filter (this is given 

by the first term in the summation on the RHD term); and penalising for poles located far from 

the unit disk, as this makes the bandwidth of filter wider and also creates bias in the resulting 

estimated frequency (this is given by the second term in the summation on the RHD term). The 

bandwidth of the notch filter is determined by p(m) as BW = vr (1 — p{m)) ifTTl . The estimated 

frequency, ^2, can be retrieved by f2 = arccos (^—^^^^ , and used to approximate the GLRT in 

(fT2l) . Although this estimator is asymptotically biased, the bias can be made arbitrarily small by 

choosing p — 1 IfTTll . However, the initial radius, p(l), should be set such that BW > 2e, to 

ensure that the realised frequency lies within the filter's range. 

Direct optimisation of (fT6l ) is difficult due to its non-linearity. However, this problem can be 
easily solved in a sequential manner, by utilising adaptive filter theory. Here we use the steepest 
descent approach to minimise the associated cost function in (fT6l) . This results in adapting (3 
and p and tracking the location of the null frequency as follows 



where p/? and pp are predefined step size parameters. Since the support of the frequency is 
bounded over the interval [ZU — e, cJ + e], the values of (3{m) can be constrained by —2 cos (cJ — e) > 
/3(m) > —2 cos (p + e). These constraints are useful in case of low SNR and prevent the filter 
from wondering outside the boundaries of the allowed frequencies. For the first sample, (m = 1), 
we center the filter frequency at /3{1) = —2 cos (uJ) and set p(l) = 1 — ^. 



Computational Complexity: here we analyse the computational complexity involved in the 
proposed algorithms. We define the overall complexity of an algorithm as the number of complex 
multiplications (CM) and complex additions (CA) required. 

Energy Detector - the test statistics is T(Yi:m) = I^m=i |y("^)l^- The computational complexity 



(3{m) = (3{m - 1) - 



p(m) = p{m- 1) - fip — 




(17b) 



(17a) 



VL Computational Complexity and Analysis 
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is therefore M (CM + CA). 

Single Matched Filter - the test statistics is defined in ( TTal ). Obtaining r (w) involves evaluating 
Uc (w) and Us (w). The total computational effort is therefore 2M (CM + CA). 
CANF based Detector - the steps involved in the proposed algorithm are: 

1) Evaluation of the filter in ([JS]): 3M (CM + CA). 

2) Adapting (5 in (flTal) : 2M (CM + CA). 

3) Adapting p in (flTbl) : 3Af (CM + CA). 

4) Performing a single matched filter: 2M (CM + CA). 

Therefore, the overall complexity of the proposed algorithm is lOM (CM + CA). 
Performance Analysis: we now provide an analysis regarding the performance gain obtained 
by the CANF algorithm over the Energy Detector and the Matched filter. 

• CANF vs. Energy Detector: here we show the SNR ratio gain of the CANF over the Energy 
detector. We do so by evaluating their ratio of SNRs 

SNRcANF Ps ,( Ps \ _ /^.Norfw _$ 



-/ 



» 1 



(18) 



SNRed Nodu \ U No du J No rfo; 2e 

where Ps is the signal's energy, Nq is the PSD of the AWGN and $ is the system's band- 
width. Thus, the narrower the notch filter is, the greater the SNR improvement, providing 
better performance. This is directly related to the second term of the optimisation formulation 
in (dH). 

Matched Filter analysis: here we demonstrate the sensitivity of the matched filter to fre- 
quency offsets. In the complex domain the test statistic, r (cJ) in (TTal) can be expressed 
as 



r (iO) \Hi 



M 



m—l 



Y('m) exp {—jmVi) 
|if||exp(je)| 



M 



(H c-xjp {jmil + G) + V{m)) cxp {—jmuj) 



■m—l 



zero mean AWGN 



M 



exp (jmn) cxp {—jmuj) 



\H\ 



S—i^ — LJ 



M 



^ exp (jmS) 



<^ 0. 

M>>2TrS 

(19) 



Hence, the test statistic converges to for frame lengths that are larger than the period of 
the frequency offset, 5, rendering the matched filter incapable of performing robust LRT. 
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VII. Simulation Results 

We present the performance of the proposed algorithm and comparison via simulations. We 
begin by evaluating the performance of the adaptive notch filter to perform frequency estimation. 
We tested the algorithm for a frame length of M = 64. The nominal angular frequency was 
set to 1.9635 and the maximal allowed offset was set to 0.98. The realised angular frequency 
was set to 2.45. Fig. [2] presents boxplot results for normalised frequency estimation error of the 
CANF, which clearly shows the good performance obtained by the notch filter based frequency 
estimator. 

Next, we compare the proposed algorithm (labeled as CANF) with other detection schemes: 
the energy detector; the "close to optimal" solution based on a bank of matched filters as per (fTTI) . 
with K = 20 (this makes the computational complexity of this algorithm roughly 4 times the 
complexity of the CANF algorithm); a mismatched detector which makes the realised frequency 
is the nominal one; and as a lower bound we use a detector that has a full knowledge of 
the realised frequency; The Receiver Operating Characteristics (ROC) results are presented in 
Fig. [3l for various SNR values. As seen from the results, the proposed algorithm performs much 
better than the energy and the mismatched detectors. The bank of matched filters provides poorer 
performance than the CANF detector. It would take K fa AO to achieve similar performance as the 
CANF detector, which makes its computational complexity around 8 times more for the same 
performance characteristics. The mismatched filter performs very poorly, which demonstrates 
how important it is to take into account frequency offsets in the design of the system. As the 
results depict, the CANF performs close to the lower bound. 

Next, we fixed PpA to 0.1 and obtained Pd for various frame lengths and different SNRs. The 
simulation results are presented in Fig. IH We observe that the mismatched detector performs 
poorly and that increasing the frame length does not improve its performance. We also note that 
the proposed algorithm performs very close to the lower bound. 

VIII. Conclusions and Future Work 

In this paper, we proposed a low complexity algorithm for spectrum sensing over fading 
channels with frequency offsets. The scheme was based on an adaptive notch filter to perform a 

February 16, 2011 DRAFT 



12 



low complexity frequency estimation followed by a single matched filter. Simulation results show 
that the comparable performance to the "close to optimal" scheme can be obtained with only a 
fraction of the algorithmic complexity. Future research will include the scenario of collaboration 
of multiple sensors, and dynamic evolution of the frequency offset. 
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Fig. 1. Spectrum of ATSC channel spectrum. The channel occupies 6 MHz and is relatively flat except for the pilot signal 
located in 310 kHz above the lower edge of the channel. 
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Fig. 2. Frequency estimation error of the adaptive notch filter Vs. SNR for A'^ = 64 
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Fig. 3. Comparison of Paid vs. Pfa for the proposed algorithm, for TV = 64 and SNR = [0, 3, 6]dB 
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Fig. 4. Detection probabilities for various frame lengths (A^ = {64, 128, 256}) and Pfa — 0.1 
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